Bulletin of Electrical Engineering and Informatics 
Vol. 9, No. 4, August 2020, pp. 1461~1467 
ISSN: 2302-9285, DOI: 10.1159 1/ee1.v914.2347 O 1461 


Multi-wavelet level comparison on compressive sensing for MRI 
image reconstruction 


Indrarini Dyah Irawati, Sugondo Hadiyoso, Yuli Sun Hariyani 
School of Applied Science, Telkom University, Indonesia 


Article Info ABSTRACT 

Article history: In this study, weproposed compressive sampling for MRI reconstruction based 
on sparse representation using multi-wavelet transformation. Comparing the 

Received Oct 24, 2019 performance of wavelet decomposition level, which are level 1, level 2, level 3, 

Revised Feb 17, 2020 and level 4. We used gaussian random process to generate measurement 

Accepted Mar 23, 2020 matrix. The algorithm used to reconstruct the image is f,norm. The 


experimental results showed that the use of wavelet multi-level can generate 
higher compression ratio but requires a longer processing time. MRI 
Keywords: reconstruction results based on the parameters of the peak signal to noise 
ratio (PSNR) and structural similarity index measure (SSIM) show that the 


Compressive sampling higher the level of decomposition in wavelets, the value of both decreases. 


Multi-wavelet 


MRI 
Reconstruction 
This is an open access article under the CC _BY-SA license. 
COO) 
Corresponding Author: 


Indrarini Dyah Irawati, 

School of Applied Science, 

Telkom University, 

Jl. Telekomunikasi, Terusan Buah Batu Bandung 40257, Indonesia. 
Email: indrarini @tass.telkomuniversity.ac.id 


1, INTRODUCTION 

Today, there are many electronic devices produce digital information such as image, voice, video 
and data. Digital imagesincluding digital photography, medical photos (magnetic resonance imaging/MRI, 
X-ray), seismic image have large size so as require big storage capacity. These images also consume a wide 
bandwidth when transmitted. MRI medical image has file size about 5-6 MB with 100-1.000 images have 
contains in series 4-19 [1]. The MRI needs a large enough capacity so that it takes a technique to reduce its 
capacity thus saving storage space but still ensure the MRI image quality. A new method introduced 
by Donoho [2] is called compressive sensing (CS) can perform the sampling and reconstruction process 
simultaneously. CS presents lower signal acquisition rate under the Nyquist rate when the signal known sparse 
without any loss quality. The signal can be expressed by few samples therefore the signal will be compressed. 

There are two properties that must be met at compressive sensing, ie sparsity [3] and incoherency [4, 5]. 
Sparsity denotes that a signal can represent in the domain with small number of non-zero samples. Incoherency 
means that there is no correlation between two different basis, ie the representation base and projection base. 
An opportunity to achieve sparse signal on the right domain so that CS performs compression with minimal 
sample with high degree of accuracy. 

Compressive sensing application usage has been done in various field such as photography image [6], 
medical image [7-10], network monitoring [11], holographic microscopy image [12], weather monitoring 
system [13], and so on. In this paper, we explore the sparsity MRI image representation using wavelet 
transformation. MRI is a research field that attract the attention of researchers until now. Some research has 
developed the CS scheme specifically designed for the MR image recovery algorithms [14, 15]. 
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Other research, study the sparse MRI representation methods [16-18]. Hot, et al. [16] compared the masked 
coefficients number of 2D discrete fourier transform (DFT) and 2D discrete cosine transform (DCT) 
as significant values of the MR image. The experiment result shows that PSNR depends on the mask which 
groups a number of high frequency coefficients. Guerquin-Kern, et al. [17] combined sampling method with 
arbitrary k-space trajectories using wavelet and iterative shrinkage/thresholding algorithm (ISTA) 
to accelerate MRI reconstruction. The proposed method is superior to total varian (TV) regularization in 
terms of image quality. In [18], authors proposed four level wavelet-encoded CS (WE-CS) to recover MR 
images using reduced samples. This method improves the performance and reduces the sampling time 
compared to fourier-encoded CS (FE-CS). MRI research on wavelet-based image reconstruction has actually 
been carried out by several researchers. Some researchers have proposed methods or algorithms for MRI 
image reconstruction specifically wavelet-based. Research by Chen and Huang [19], proposed a wavelet tree 
structure to improve the performance of conventional compressed sensing MRI. Research by Ragab et. al [20], 
proposed compressive sensing MRI reconstruction using empirical wavelet transform. This research claims 
that the proposed method outperforms the state-of-the-art methods in terms of signal-to-noise ratio and 
structure similarity metrics. Similar research by Ragab et. al [21], simulates the dual tree complex wavelet 
transform method on CS-MRI. This research has a better SNR performance compared to MRI wavelet tree 
sparsity. Lai et. al [22] has proposed compressed sensing MRI using graph-based redundant wavelet 
transforms. However, this proposed method is only able to produce an SNR of less than 20 dB at a 20% 
sampling rate. Zhu et. al [23] has proposed CS MRI reconstruction using a DWT dual-density dual-tree 
complex. This method modifies conventional wavelet decomposition to get more detailed sub-bands that are 
expected to reduce aliasing artifacts. These reviews provide a thought that CS techniques will become an 
interesting topic in the future to overcome the limitations of memory storage. The main issue and also the 
problem of implementing CS especially on images is how to get a high compression ratio and PSNR with 
simple techniques and produce high accuracy. This paper proposes a compressive sensing with 
the multi-wavelet technique which is applied to MRI images. This method is combined with the ?;norm to 
reconstruct the image algorithm. This proposed method becomes the state of the art of the proposed method 
in related research above. The aim is to analyze the effect of multi wavelets on the performance of the 
proposed method so that the best approach can be applied in the future. The structure of the paper is arranged 
as follows, section 2 presents the proposed method for MRI reconstruction system. The experimental results 
and analysis are provided in section 3. In section 4 is shown summary and conclusions, also further research. 


2. RESEARCH METHOD 

The proposed of compressive sensing simulation modelingfor the MRI reconstruction is shown in 
Figure 1. MRI input is an image with format extensions such as .png and .bmp. Sparsity base aims to produce 
Sparse image which has non-zero values mouth slightly. In this research, we explore the sparsity by 
comparing the level of wavelet transformation. We also examine the effect of level wavelet transformation 
such as level 1, level 2, level 3 and level 4. We take a M X N mesurement matrix (®) that produced from 
random gaussian process. Reconstruction algorithms serves to restore the compressed image. We use 
£,norm for reconstructing [24]. MRI output is the result of the inverse wavelet transform. 

In DWT 2D image used by the filter process using wavelet bases in rows and columns that make 
the matrix approximation and detail. The decomposition process is tiered so-called multi wavelet as shown in 
Figure 2. In the matrix approximation done back in the filter process to produce a matrix approximation and 
detail at the next level. In this study a comparison effect level wavelet used against the parameters measured. 
Wavelet is a more efficient method for estimating original images with a small number of basic elements. 
In mathematical equations, wavelet transform on the function x (t) is interpretation of these functions into 
a set of basis functions below. 





Sparsity Base (Y): 
(Exploration of Measurement Reconstruction Inverse 
—+| Mother Wavelet Base (®): | Algorithm: Wavelet 
MRI and Level Wavelet random gaussian £,norm preustocm MRI 
Input Transformation) ie i, Output 





Figure 1. The proposed of compressive sensing simulation modelingfor the MRI reconstruction 


X(t) = JX, x(t). Yap (t) at (1) 
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where Y, »(t) is a basic function, which is a bandpass version of an enlarged and shifted time signal '¥ (t). 
The signal’? (t)is named the mother wavelet and expressed as follows: 


Van = = ¥(—) (2) 


a 


where a is the enlargement parameter and b is the shift parameter. 
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Figure 2. 2D wavelet decomposition [25] 


Donoho introduces compressive sensing as a new method of signal processing [2]. The main idea 
is to take the signal with less sample below Nyquist theory where the sampling results meet the original 
sparse signals. The original has a sparse representation thereof if it 1s expressed in the appropriate base. 
The mathematical equations for a vector xe€R”is represented using the wavelet basis V = [W,W¥.... Py], 
vector x can be formulated [11]. 


f= sO) 3) 


wheres is the coefficient of x obtained from s; = (x,¥,). 
In matrix form can be written: 


x= Ws (4) 


where Y an N x N matrix, x and s is acolumn vector N x 1. 

Signals x expressed K-sparse if only K from coefficients of s, when the value is not zero while 
the (N-K) coefficient is zero. In compressive sensing is required M sample of the signal, where K <M«N. 
As shown in Figure 1, y is a measurement vector with length (M), then can be expressed in the equation: 


y= ®x=OVs (5) 


with y is the output signal that has been compressed, ® is the projection matrix with size M x N, ¥ , a base 
matrix with size N x N, and 9 = ®wW 1s the matrix Mx N. 

Based on the equation 5, to get the compressed signal y, the required measurement matrixM x N. 
The measurement matrix used is a gaussian random matrix. The smaller the value of M (measurement rate) then 
the resulting matrix size is also smaller. In this study measured the effect of the value of M or the number of 
sampling of the measured parameters namely compression ratio (CR). 

After getting the uncompressed signal (y), then do the reconstruction process to obtain the original 
image. Some reconstruction algorithm that can be used include norm-order one, orthogonal matching pursuit 
based greedy algorithm, and others.In this experiment, the algorithm first order norm for reconstruction 
poses. Image reconstruction results will be compared with the original image and then calculated PSNR and 
SSIM. Another measurement parameter is the duration of the process named processing time. 
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3. RESULTS AND ANALYSIS 

In this study, we measured several parameters, such as PSNR, SSIM, CR, and processingtime 
for the use of wavelet and multi wavelet. Matrix measurement used is a gaussian random matrix with testing for 
M different values of 10% to 90%. While the reconstruction algorithm used is £;norm. Tests performed on the 
images consist of cephal MRI with a size of 256 x 256 and brain MRI transversal with size 350 x 227. Because 
the matrix is used random matrix, then tested 5 times and then calculated the average for each parameter. 

The PSNR parameter test results can be seen in Figure 3. The test is performed for cephal MRI. 
PSNR illustrated the quality of the reconstructed image. The test results show that the greater the measurement 
rate,the PSNR value increases. Experiments at the DWT decomposition level show that the more decomposition 
levels, the PSNR is decreasing. According to [26], the image quality is declared good when 
it has a PSNR= 30dB. According to Figure 3, DWT level 1 and level 2 have good image quality 
at measurement rate > 20%, whereas DWT level 3 and level 4 have good quality at measurement rate > 80%. 

SSIM is a parameter that shows the level of resemblance between the original image and 
the reconstructed image. Figure 4 presents the SSIM values at different DWT decomposition levels on cephal 
MRI, which illustrates the greater the value of the measurement rate the higher the SSIM value. From the test 
results, level 1 decomposition has the highest SSIM compared to level 2, level 3, and level 4, almost all 
of the measurement rate. 

CR is used to show the performance of the compression algorithm in the relative reduction in the size 
of the resulting data compressionwhich expressed as uncompressed size division by compressed size. 
Figure 5 shows CR on cephal MRI, which illustrates that the greater the measurement rate, the smaller the CR. 
When compared between levels of DWT decomposition, the greater the level of decomposition, the CR is 
increasing. In this case, CS on DWT with a high level of decomposition is suitable for obtaining large CRs thus 
saving storage space. 

Figure 6 shows the processing time for cephal MRI. The processing time is proportional 
to the complexity of the algorithm 1s executed. Multi tiered process on wavelet scheme requires a longer time 
than the wavelet level 1. The processing time at DWT level 4 is the longest compared to the others. 

Based on PSNR, SSIM, CR, and processing time metrics it can be explained that PSNR is proportional to 
SSIM and inversely proportional to CR and processing time. The higher level of DWT decomposition will result in 
high CR with the risk of a long processing time, but reduces the accuracy seen from the PSNR and SSIM 
parameters. In these experiments, the DWT decomposition level 2 is suitable for obtaining high CR while still 
obtaining good reconstruction image quality with a minimum measurement rate of 20%. 
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Figure 5. CR for measurement rate [10% -90%] Figure 6. Time of processing for measurement rate 

at different DWT decomposition levels [10% -90%] at different DWT decomposition levels 


Figure 7 and Figure 8 show the comparison between the original image with the reconstruction image 
at DWTdecomposition level 1 (b) and (c) and DWT decomposition level 2 (d) and (e) for measurement rates 
of 50% and 90%. Figure 7 for cephal MRI while Figure 8 for brain MRI transversal. According to Table 1, 
for measurement rate 50% in cephal MRI with DWT decomposition level 2 is obtained PSNR=38.12dB, 
SSIM=0.486, and CR=1.85. When seen in plain view, the results of reconstruction at measurement rate=50% 
are still unclear, different from the results of DWT decomposition level 1 for the same measurement rate. 
Image obscurity was also seen in experiments on brain MRI transversal, with the value of each parameter 
as follows: PSNR= 39.89dB, SSIM=0.499 and CR=?2.11. 





(d) (e) 


Figure 7. Cephal MRI 256x256, (a) Original image, (b) Reconstruction image using level | (measurement 
rate=50%), (c) Reconstruction image using level 2 (measurement rate=90%), (d) Reconstruction image using 
level 2 (measurement rate=50%), (e) Reconstruction image using level 2 (measurement rate=90%) 


Multi-wavelet level comparison on compressive sensing for MRI image... (Indrarini Dyah Irawati) 


1466 O ISSN:2302-9285 


4. 





(d) (e) 


Figure 8. Brain MRI transversal 350x227, (a) Original image, (b) Reconstruction image using level 1 
(measurement rate=50%), (c) Reconstruction image using level 1 (measurement rate=90%), 


(d) Reconstruction image using level 2 (measurement rate=50%), (e) Reconstruction image using level 2 


(measurement rate=90%) 


Table 1. PSNR, SSIM, CR for different level decomposition for measurement rate=50% 


aes Cephal MRI Brain MRI 
PSNR (dB) SSIM = CR  PSNR(dB) SSIM CR 
Level 1 44.59 0.689 2.76 47.76 0.712 3.01 
Level 2 38.12 0.486 1.85 39.89 0.499 2.11 
Level 3 28.19 0.243 1.46 32.97 0.313 1.54 
Level 4 28.15 0.198 1.23 32.67 0.245 1.33 


CONCLUSION 
In this paper, several conclusions are, the higher the number of measurement rate the greater PSNR 


and SSIM value are generated, the smaller the CR is obtained, and the longer the processing time is needed. 
The use of multi wavelet scheme is able to increase the CR but reduce PSNR and SSIM and delayed 
processing time. In this experiment, DWT level 2 is recommended for compression with a measurement rate 
> 50% to get a good MRI reconstruction quality. 


Future research is to develop a method for obtaining sparsity based on wavelet technique by forming 


a new matrix containing wavelet coefficients with significant values. In addition, we will explore the effects 
of using other measurement matrices such as Bernoulli, Binomial, the use of proposed methods for color 
MRI or dynamic MRI, observing the effect of using other reconstruction algorithms. 
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